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EXPERIMENTAL  RANDOM  FATIGUE  IN  ELASTIC-PLASTIC  RANGE- 
FIRST  ORDER  MODELS 

T.  C.  Huang  and  Vinod  K.  Nagpal 
Department  of  Engineering  Mechanics 
University  of  Wisconsin-Madison 
Madison,  Wisconsin  53706 

ABSTRACT 

An  experiment  based  on  11  probabilistic  parameters  and  experiment 
design  is  conducted  for  fatigue  of  materials  in  an  elastic-plastic  range 
under  random  vibrations  to  develop  regression  models  which  will  predict 
the  fatigue  life.  The  experiment  consisted  of  a  total  of  24  tests  and 
is  conducted  in  3  designs.  The  response  of  each  test  specimen  is  digi¬ 
tized,  recorded,  and  is  fitted  with  time  series  models.  The  time-series 
model  parameters  are  then  used  to  compute  the  spectrum  and  spectral 
moments  from  which  the  magnitude  levels  of  probabilistic  parameters  are 
computed  and  coded.  The  first -order  regression  models  are  developed 
for  each  design  by  regressing  the  log  of  fatigue  life  on  the  corre¬ 
sponding  coded  parameters.  The  tables  of  the  analysis  of  variance,  and 
of  the  predicted  lives  together  with  residuals  and  952  confidence  inter¬ 
vals  are  constructed  for  each  model.  From  the  analysis  of  variance  the 
F-ratio  is  computed  to  check  whether  the  model  is  acceptable.  The  first- 
order  model  Involving  all  11  parameters  based  on  24  tests  is  considered 
to  be  the  statistically  best  one.  The  methodology  of  this  report  is 
found  to  be  an  accurate  and  reliable  approach  and  is  suitable  for 
industrial  applications.  This  method  constrasts  with  the  prevailing 
method  which  is  based  on  linear  damage  accumulation  and  cycle  counting, 
involves  a  several  hundred-percent  error  as  a  rule,  and  is  hardly 
applicable  to  industrial  problems. 
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INTRODUCTION 

The  principle  of  linear  damage  accumulation  based  on  cycle  counting 
has  been  consistently  found  to  give  unreliable  and  inaccurate  estimates 
of  fatigue  life  under  random  vibrations  [1-10].  This  group  of  represen¬ 
tative  references  has  been  discussed  in  [11]. 

Recently  a  new  history-dependent,  stochastic  model  of  cumulative 
damage  is  being  developed  by  Bogdanoff  [12-14]  by  taking  a  comprehen¬ 
sive  view  of  the  entire  failure  process  to  improve  the  predictive 
accuracy.  Another  approach,  entirely  different  from  all  above  mentioned 
approaches,  has  been  developed  in  [11,15]  to  predict  the  fatigue  life 
of  materials  under  random  vibrations  in  the  elastic  range.  The  same 
approach  is  used  now  for  elastic-plastic  range.  An  experiment  program 
of  24  tests  based  on  11  probabilistic  parameters  and  experiment  design 
is  conducted  in  3  designs.  For  each  design  the  regression  models  are 
developed  to  predict  fatigue  life  under  random  vibrations.  Again,  this 
approach  has  been  found  to  give  very  reliable  and  accurate  estimates. 

The  purpose  of  this  research  is  to  develop  a  novel  methodology 
for  random  fatigue  based  on  the  probabilistic  parameters  and  experiment 
design.  Therefore,  the  choice  of  the  specimen  materials,  types  of 
loadings  (axial,  shear  or  bending),  and  range  of  parameters  is  irrele¬ 
vant  and  they  are  arbitrarily  chosen.  p — 

.  EXPERIMENTS,  PARAMETERS.  DESIGNS  AND  MODELS 


1.  Experiment 

The  experiment  program  was  conducted  with  24  specimens  in  3 
designs.  These  specimens  were  numbered  from  25-48  to  distinguish 
them  from  specimens  numbered  from  1-24  used  in  the  elastic  case  [11].  d-?8 

)Dj-t  j  special 
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The  bending  specimen  of  aluminum  alloy  6061-T6  and  Che  experimental  set¬ 
up  described  in  [11,16,17]  were  also  used  for  this  experiment. 

2.  Parameters 

A  typical  random  response  signal  in  elastic-plastic  range  which 
is  characterized  by  11  probabilistic  parameters  is  given  in  Fig.  1. 

They  are  mean,  variance,  zero  upcrossings,  level  upcrossings,  Cy 
level  upcrossings,  duration  of  excursion  above  zero,  duration  of  excur¬ 
sion  above  £,  level,  duration  of  excursion  above  e  level,  band 
f  y 

width,  average  amplitude  above  level,  and  average  amplitude  above 

£y  level.  In  these  parameters  £^  represents  the  strain  level  corre¬ 
sponding  to  a  material  life  of  10^  cycles  in  the  deterministic  fatigue 
tests,  and  £y  represents  the  yield  point  strain  of  the  material.  A 
mathematical  analysis  of  these  11  probabilistic  parameters  is  given  in 
[18].  Hereafter,  these  probabilistic  parameters  will  be  simply 
referred  to  as  variables  as  well. 

3.  Designs 

The  first  design  of  the  experiment  consisted  of  10  tests  which 
formed  a  full  factorial  design  with  2  center  points.  The  second  design, 
a  central  composite  design  with  4  center  points, was  formed  by  adding  8 
more  tests  to  the  first  design.  The  third  design  was  formed  by  adding 
to  the  second  design  random  replications  of  six  tests  from  the  second 
design. 

a.  Range  The  first  design  was  based  on  the  three  different 


magnitude  levels  of  the  3  variables,  namely,  mean,  variance  and  zero 
upcrossings.  Two  more  levels  from  each  of  the  3  variables  were  added 
in  the  second  design  so  that  one  was  lower  and  the  other  one  higher 
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Chan  the  3  previous  levels.  In  Che  chlrd  design,  6  tests  were  repli¬ 
cated  without  introducing  new  levels  of  the  variables.  Thus,  the 
second  and  the  third  designs  had  5  levels  from  each  of  the  above  3 
variables.  The  lowest  and  the  highest  of  the  5  levels  represent  the 
lower  and  upper  bounds  of  the  range  of  variables.  The  upper  bounds  of 
the  mean  and  variance  are  limited  such  that  the  maximum  capability  of 
the  shaker  table  displacement  is  not  exceeded.  Whenever  the  displace¬ 
ment  exceeds  the  maximum  capability  of  the  shaker,  the  overtravel 
limit  is  reached  and  the  shaker  is  shut  off  automatically.  Otherwise, 
the  range  of  the  variables  can  be  arbitrarily  chosen.  The  ranges  of 

mean,  variance  and  zero  upcrossings  based  on  5  levels  are  241-2657 

6  6  2 
micro  inches/inch,  2.42x10  -  4.36x10  (micro  inches/inch)  ,  and 

9.12  -  15.32,  respectively. 

b.  Code  The  above  3  variables  which  were  controlled  in  the 
experiment  were  coded.  The  remaining  8  variables  which  could  not  be 
controlled  in  the  experiment  were  computed  from  the  random  response  of 
each  test  specimen  as  described  in  [11].  The  variables  and  the  actual 
levels,  and  the  coded  variables  and  coded  levels  are  shown  in  Table  1. 
Tables  2  and  3  show  the  actual  levels  of  variables  and  coded  levels  of 
coded  variables,  respectively,  of  all  11  variables  for  the  24  tests. 

4.  Models 

The  regression  models,  also  referred  to  as  life  predicting  equa¬ 
tions,  were  obtained  for  each  of  the  3  designs  as  was  donp  in  [11].  For 
the  first  design  only  one  life  predicting  equation  of  3  variables  was 
obtained.  For  each  of  the  second  and  third  designs  two  life  predicting 
equations  were  obtained,  one  Involves  only  3  variables  and  the  other 


one  involves  all  11  variables.  These  equations  were  obtained  by 
regressing  the  log  of  actual  life,  y,  on  the  corresponding  coded  vari¬ 
ables  involved  in  the  model.  The  models  predicted  the  log  of  life,  y. 
The  predicted  fatigue  life,  T,  is  computed  by  talcing  the  antilog  of  y. 

The  table  of  analysis  of  variance,  and  the  table  of  predicted 
lives  together  with  residuals  and  95%  confidence  intervals,  are  con¬ 
structed  for  each  model.  From  the  analysis  of  variance,  the  F-ratlo 
for  each  model  was  computed  to  determine  whether  the  model  Is  accept¬ 
able.  The  confidence  intervals  were  computed  using  the  standard 
deviations  of  y  and  the  t  values  from  the  t-table  with  degrees  of 
freedom  equal  to  that  of  residuals.  The  distribution  of  residuals  was 
studied  for  any  pattern  or  trend  present.  If  the  pattern  or  trend 
existed,  the  model  was  considered  to  have  missed  some  significant 
variables.  The  analysis  of  3  designs  is  given  below. 

U.  FULL  FACTORIAL  DESIGN  WITH  TWO  CENTER  POINTS 

The  first  ten  tests  from  25-34  of  the  24  tests  formed  a  full 
factorial  design  with  2  center  points.  One  life  predicting  equation 
of  three  variables  is  obtained  by  regressing  log  of  actual  life,  y, 
on  coded  levels  of  mean,  variance  and  zero  upcrossings.  These  three 
variables  are  represented  by  symbols  x^-x^.  The  life  predicting 
equation  of  all  11  variables  could  not  be  obtained. 

1.  Three  Variables 

The  life  predicting  equation  is  obtained  as 

9  -  3.52  +  0.0538x3^  -  0.341X2  -  0.0514x3  (1) 

The  analysis  of  variance  of  equation  (1)  is  given  in  Table  4.  The  F- 
ratlo  of  4.328  with  3  and  6  degrees  of  freedom  is  found.  The 
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corresponding  F  value  from  the  F- table  is  4.76  which  is  greater 
than  the  F-ratio  of  the  equation.  This  inplies  that  the  regression 
is  not  effective  and  the  model  is  not  acceptable.  The  residual  sum  of 
squares  is  0.4862  in  comparison  to  a  total  of  2.0269,  a  32.6%.  The 
other  68.4%  of  the  total  sum  of  squares  is  due  to  regression. 

This  model  does  not  qualify  the  F-test  and  may  not  be  acceptable. 
Nevertheless  the  predicted  lives  together  with  their  residuals  and  95% 
confidence  intervals  are  computed  as  given  in  Table  5  in  order  to  see 
how  the  results  of  an  unacceptable  model  appear.  The  confidence 
intervals  are  fairly  wide  because  the  t  value  with  six  degrees  of 
freedom  is  high.  The  actual  life  of  the  test  number  30  only  is  slightly 
below  the  lower  limit  of  the  predicted  interval. 

2.  All  11  Variables 

The  life  predicting  equation  based  on  all  11  variables  could  not 
be  obtained  because  the  number  of  tests  in  this  design  is  smaller  than 
the  variables  involved. 

m .  CENTRAL  COMPOSITE  DESIGN  WITH  FOUR  CENTER  POINTS 

The  full  factorial  design  is  augmented  by  8  more  tests  to  form  a 
central  composite  design  with  4  center  points  resulting  in  a  total  of 
18  tests.  Two  life  predicting  equations  are  obtained  for  this  design. 
The  first  one  consists  of  the  same  3  variables,  while  the  second  one 
consists  of  all  11  variables. 

1.  Three  Variables 

The  life  predicting  equation  is  obtained  as 
9  -  3.55  +  0.0242x1  -  0.327x2  -  0.0597x3 


(2) 


7 


The  analysis  of  variance  of  equation  (2)  is  given  in  Table  6.  The  F- 
ratio  of  14.80  with  3  of  14  degrees  of  freedom  is  obtained.  The  corre¬ 
sponding  F  value  from  the  F-table  at  95%  significance  level  is  3.34. 
The  table  value  of  F  is  smaller  than  the  F-ratio  which  implies  that 
the  regression  is  effective  and  the  model  is  acceptable.  The  residual 
sum  of  squares  is  0.4862  in  comparison  to  a  total  of  2.0269,  a  24.0%. 

The  other  76.0%  of  the  total  sum  of  squares  is  due  to  regression.  In 
this  case  the  sum  of  squares  due  zero  upcrossings  is  higher  than  the 
sum  of  squares,  due  to  mean. 

The  predicted  lives  together  with  residuals  and  95%  confidence 
interval,  are  given  in  Table  7.  The  confidence  intervals  are  fairly 
narrow.  The  actual  lives  of  lesL  numbers  30,  33  nnd  34  are  slightly 
off  from  the  corresponding  predicted  interval.  The  residuals  also  an 
relatively  large  in  magnitude.  The  investigation  of  this  model  suggests 
that  more  variables  which  have  significant  effects  on  the  fatigue  life 
should  be  introduced. 

2.  All  11  Variables 

The  life  predicting  equation  is  obtained  as 

y  *  3.39  -  0.0143x^  -  0.463x2  -  0.144x^  +  0.020x^ 

-  0.314x^  -  0.727x^  +  0.637x2  +  0.112xg  -  0.0103x^ 

+  0.101x10  -  0.0803x11  (3) 

The  analysis  of  variance  of  equation  (3)  is  given  in  Table  8.  The  F- 
ratio  of  6.036  with  11  and  6  degrees  of  freedom  is  obtained  for  this 
model.  The  corresponding  F  value  from  the  F-table  is  2.92  at  95% 
significance  level.  The  F-ratio  is  greater  than  the  table  value  of  F 


which  means  that  the  regression  is  effective  and  that  the  model  is 


acceptable.  The  sum  of  squares  due  to  residuals  is  0.1682  in  compari- 

sion  to  a  total  of  2.0269,  a  8.3%.  The  other  91.7%  of  the  total  sum 

* 

of  squares  is  due  to  regression.  The  sum  of  squares  due  to  6  of  the 
11  variables  is  considered  insignificant.  These  six  variables  are 


mean,  level  upcrossings,  duration  of  excursion  above  zero  and  l 


levels,  band  width  and  average  amplitude  above  level. 


The  predicted  lives  together  with  residuals  and  95%  confidence 
Intervals  arc  given  in  Table  9.  The  confidence  intervals  are  relatively 
wide  because  t  value  with  6  degrees  of  freedom  associated  with  the 
residuals  is  high.  The  actual  lives  of  all  tests  fall  within  the 
predicted  confidence  intervals.  The  residuals  are  small  in  magnitudes. 

IV.  CENTRAL  COMPOSITE  DESIGN  WITH  FOUR  CENTER  POINTS  AND  SIX  REPLICATIONS 


Six  more  tests  are  added  to  the  central  composite  design  by  randomly 
replicating  six  of  the  18  tests.  Adding  these  six  tests  resulted  in 
a  total  of  24  tests  in  the  design.  Two  life  predicting  equations  are 
obtained  for  this  design  also,  one  using  1  variables  and  another  one 
using  all  11  variables.  The  equal  ions  and  their  analyses  are  as  ioilnw.s. 
1.  Three  Variables 

The  life  predicting  equation  is  obtained  as 

y  -  3.5C  i  0.0195x1  -  0.125*.,  -  O.UilSx^  (4) 

The  analysis  of  variance  oi  cquai  *  ■  -r  in  Table  iO.  The  r- 

ratio  for  the  equation  above  '•>  i.  .  and  20 

degrees  of  freedom  which  is  irea.fi  ’tun  ttu  to  -...ui.i,  .  value 

of  3.10  from  the  K-tubi«  at  95  X  *  . .  ti.ar.  <  ...  i  i<  »  that 

the  regression  is  effective  and  ttu-  m»>  cl  is  .«t  .e,.r  .<o  i«.  r,e  residual 
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sum  of  squares  is  0.6192  in  comparison  to  a  total  sum  of  squares  of  3.1801, 
a  20.5%.  The  other  79.5%  of  the  total  is  due  to  regression. 

The  predicted  lives  together  with  residuals  and  95%  confidence 
intervals  are  given  in  Table  11.  The  sum  of  squares  due  to  zero 
upcrossings  in  this  case  is  lower  than  that  due  to  mean.  The  confidence 
intervals  are  comparatively  narrow  but  the  actual  lives  of  test  nos.  27, 

29,  30,  33,  34,  35,  44  and  45  fall  out  of  the  corresponding  predicted 
confidence  intervals.  The  residuals  are  relatively  large  which  imples 
that  some  variables  which  may  have  significant  effects  on  the  fatigue 
life  are  missing  from  the  model.  The  life  predicting  equation  using 
more  variables  is  discussed  below. 

2.  All  11  Variables 

The  life  predicting  equation  is  obtained  as 

y  =  3.50  -  0.0105x1  -  0.398x2  -  0.0748x3  +  0.0459x4 

-  0.111x5  -  0.317xfi  +  0.290x?  +  0.0091xg  -  0.0383x9 
+  0.0739x10  +  0.0231xn  (5) 

The  analysis  of  variance  based  on  24  tests  of  equation  (5)  is  given 
in  Table  12.  The  F-ratio  for  the  above  equation  is  computed  to  be 
8.812  with  11  and  12  degrees  of  freedom  which  is  larger  than  the  corre¬ 
sponding  F  value  of  2.73  from  the  F-table  at  95%  significance  level. 

The  comparison  of  two  F  values  indicates  that  the  regression  is 
effective  and  the  model  is  acceptable.  The  residual  sum  of  squares 
is  0.3498  in  comparison  to  a  total  sum  of  squares  of  13.1801,  a  11.0%. 

The  other  89%  of  the  total  is  due  to  regression,  which  is  significantly 
high.  The  six  of  the  11  variables  are  considered  to  be  insignificant 
because  they  contribute  a  very  low  sum  of  squares  to  regression.  These 
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6  variables  are  level  upcrossings,  duration  of  excursion  above 

zero  level,  duration  of  excursion  above  level,  band  width,  and 

average  amplitudes  above  and  levels. 

The  predicted  lives  together  with  residuals  and  95%  confidence 
intervals  are  given  in  Table  13.  The  confidence  intervals  are  fairly 
narrow.  Actual  lives  of  all  the  24  tests  fall  within  the  predicted 
confidence  intervals.  The  residuals  appear  to  be  relatively  small  and 
their  plot  is  given  in  Fig.  2.  The  plot  shows  that  the  residuals  are 
concentrated  more  on  the  positive  side  with  fairly  small  magnitudes. 

V.  DISCUSSIONS  AND  CONCLUSIONS 


On  the  basis  of  the  analysis  of  all  models  of  equations  from  (1) 
to  (6)  it  is  observed  that  there  are  5  variables  which  have  significant 
effects  on  the  fatigue  life.  These  5  significant  variables  are  mean, 
variance,  zero  upcrossings,  level  upcrossings,  and  the  duration  of 
excursion  above  level.  The  remaining  6  variables  are  considered 

to  be  less  significant. 

Among  the  5  significant  variables  the  variance  has  the  most  sign!: i- 
cant  effect.  This  is  obvious  because  the  variance  is  a  measure  of 
square  of  the  amplitude  of  the  strain  signal.  Consequently  a  large 
variance  of  the  signal  means  large  amplitude  of  the  strain.  The  magni¬ 
tudes  of  effects  of  the  remaining  four  significant  variables  are  low 
as  compared  to  that  of  the  variance  but  are  not  negligible.  Variance, 
level  upcrossings,  and  duration  of  excursion  above  level  show 

significant  effects  in  all  3  analyses  of  10,  18  and  24  tests.  The  effect 
of  zero  upcrossings  is  found  to  be  pronounced  in  the  analysis  of  10  tests 
and  the  effect  of  mean,  in  the  analysis  of  24  tests.  Nevertheless 
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both  mean  ami  zero  tiperossings  have  been  Included  in  the  list  al  'j 
significant  variables. 

A  comparison  of  all  the  first  order  models  investigated  in  this 
report  with  respect  to  their  percent  deviations  of  the  predicted  lives 
from  the  actual  ones,  and  residual  sum  of  squares  as  percent  of  the 
total  is  given  in  Table  14.  This  table  shows  that  the  models  of  all 
11  variables  for  both  18  and  24  tests  give  a  lower  percent  deviations 
and  a  lower  percent  residual  sum  of  squares  in  comparison  to  the  model 
of  3  variables  for  the  same  design.  Based  on  24  tests,  the  model  of 
equation  (5)  of  all  11  variables  is  considered  the  best  first  order 
model  to  predict  fatigue  life  under  random  vibrations  in  the  elastic- 
plastic  range.  The  percent  deviations  of  the  best  model  range  from 
-26.0%  to  23.2%  with  an  average  deviations  of  13.4%  on  the  negative 
side  and  8.4%  on  the  positive  side.  These  deviations  are  negligible 
as  compared  to  several  hundred  percent  deviations  obtained  in  linear 
damage  accumulation  theory. 

The  lives  predicted  by  the  best  first  order  model  are  expected 
to  be  the  mean  lives  of  the  material  for  the  given  levels  of  variable.. 
The  95%  confidence  intervals  give  the  range  which  is  expected  to 
include  the  actual  fatigue  life  for  95%  of  the  tests.  The  actual  life 
of  almost  all  the  tests  is  greater  than  the  lower  level  of  confidence 
interval  so  this  level  can  be  considered  a  conservative  estimate  of 
fatigue  life. 

A  few  first  order  and  second  order  models  of  significant  variables 


will  be  given  in  [19]. 


VI.  SUMMARY 


(1)  An  experiment  of  24  tests  based  on  11  variables  and  experi¬ 
ment  design  has  been  performed  to  predict  fatigue  life  under  random 
vibrations. 

(2)  The  experiment  was  conducted  in  3  designs.  For  the  first 
design  one  first  order  life  predicting  equation  of  3  variables  has  been 
developed.  For  each  of  the  other  two  designs  two  first  order  life 
predicting  equations  of  3  and  all  11  variables,  respectively,  have 
been  developed. 

(3)  The  tables  of  analysis  of  variance,  and  of  the  predicted 
lives  together  with  residuals  and  95%  confidence  intervals  are  con¬ 
structed  for  each  equation. 

(4)  From  the  analysis  of  variance  the  F-ratio  is  computed  to 
check  whether  the  regression  is  effective  and  the  model  is  acceptable. 

(5)  Five  out  of  11  variables  have  been  found  to  have  significant 

effects  on  the  fatigue  life.  These  5  variables  are  mean,  variance, 
zero  upcrossings,  level  upcrossings  and  duration  of  excursion 
above  level. 

(6)  The  first  order  life  predicting  equation  of  all  11  variables 
and  24  tests  was  found  to  be  the  statistically  best  one. 

(7)  The  percent  deviations  of  the  predicted  lives  obtained  by 
the  best  equation  are  within  a  range  from  -26.0%  to  23.2%  with  an 
average  of  13.4%  on  the  negative  side  and  8.4%  on  the  positive  side. 
These  results  are  in  contrast  with  those  of  linear  damage  accumulation 
principle. 
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Table  1  Actual  Levels,  Coded  Levels  and  Coded  Variables 


No. 

Actual  Level 

Coded  Level 
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variaoxes 
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Center 

High 
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800.00 
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0 
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-1 

0 

+1 
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-1 

0 

+1 
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-1 
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-1 

0 
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Level  y 
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-1 
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Table  4 


Analysis  of  Variance  of  10  Testa 
First  Order  Model  of  3  Variables 
Life  predicting  equation: 


y  =  3.52 

+  0.0538^ 

-  0.341x2  - 

0.0514x3 

Source 

Sum  of 
Squares 

Degrees  of 
Freedom 

Mean 

Square 

F-Ratio 

Due  to  Mean 

0.0100 

1 

0.0100 

Due  to  Variance 

0.8439 

1 

0.8439 

Due  to  Zero  Upcrossings 

0.0264 

1 

0.0264 

Due  to  Regression 

0.8857 

3 

0.2952 

Residuals 

0.4094 

6 

0.0682 

4.328 

Total 

1.2951 

9 

F-ratio  is  smaller  than  the  table  value  of  4.76  with  3  and  6 
degrees  of  freedom  at  95%  significance  level.  So  regression  is  not 
effective  and  the.  model  is  rejected. 
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Table  6.  Analysis  of  Variance  of  18  Tests 
First  Order  Model  of  3  Variables 
Life  predicting  equation: 


y  -  3.55 

+  0.0242x1 

-  0.327x2  - 

0.0597x3 

Source 

Sum  of 
Squares 

Degrees  of 
Freedom 

Mean 

Square 

F-Ratio 

Due  to  Mean 

0.0108 

1 

0.0108 

Due  to  Variance 

1.4728 

1 

1.4728 

Due  to  Zero  Upcrossings 

0.0572 

1 

0.0572 

Due  to  Regression 

1.5407 

3 

0.5136 

Residuals 

0.4862 

14 

0.0347 

14.80 

Total 

2.0269 

17 

F-ratio  is  greater  than  the  F  value  of  3.34  from  the  table  with 
3  and  14  degrees  of  freedom  at  95%  significance  level.  So  regression 
is  effective  and  the  model  is  accepted. 


Table  7.  Results  of  18  Tests,  First  Order  Model  of  3  Variables 
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Table  8.  Analysis  of  Variance  of  18  Tests 
First  Order  Model  of  11  Variables 

Life  predicting  equation: 


y  =  3.39  -  Q.0143x^  -  0.463x2  -  0.144x2  -  0.0200x^ 

-  0.314x2  -  0.727Xg  +  0.637x2  +  0.112xg  +  0.0103xg 
+  0.101x1Q  -  0.0803xn 


Source 

Sum  of 
Squares 

Degrees  of 
Freedom 

Mean 

Squares 

F-Ratio 

Due  to  Mean 

0.0108 

1 

0.0108 

Due  to  Variance 

1.4728 

1 

1.4728 

Due  to  Zero  Upcrossings 

0.0572 

1 

0.0572 

Due  to  Level  Upcrossings 

of  Level 

0.0009 

1 

0.0009 

of  e  Level 

0.0794 

1 

0.0794 

Due  to  Duration  of  Excur- 

sion  Above  Zero  Level 

0.0349 

1 

0.0349 

Due  to  Duration  of  Excur- 

sion  Above  Level 

0.1650 

1 

0.1650 

Above  E  Level 

y 

0.0015 

1 

0.0015 

Due  to  Band  Width 

0.0002 

1 

0.0002 

Due  to  Average  Amplitude 

Above  ef  Level 

0.0097 

1 

0.0097 

Above  £  Level 

y 

0.0262 

1 

0.0262 

Due  to  Regression 

1.8587 

11 

0.1690 

Residuals 

Total 

0.1682 

2.0269 

6 

17 

0.0280 

6.036 

F-ratio  is  greater  than  the  table  value  of  2.92  with  11  and  6  degrees 
of  freedom  at  95%  significance  level.  So  regression  is  effective, 
and  the  model  is  accepted. 


Table  10.  Analysis  of  Variance  of  24  Tests 
First  Order  Model  of  3  Variables 
Life  predicting  equation: 


y  =  3.56 

+  O.OlOSx^^ 

-  0.325x2  - 

0.0518x3 

Source 

Sum  of 
Squares 

Degrees  of 
Freedom 

Mean 

Squares 

F-Ratio 

Due  to  Mean 

0.1223 

1 

0.1223 

Due  to  Variance 

2.3905 

1 

2.3905 

Due  to  Zero  Upcrossings 

0.0482 

1 

0.0482 

Due  to  Regression 

2.5609 

3 

0.8536 

Residuals 

0.6192 

20 

0.0310 

27.54 

Total 

3.1801 

23 

F-ratio  is  greater  than  the  table  value  of  3.10  with  3  and  20  degrees 
of  freedom  at  95%  significance  level.  So  regression  is  effective,  and 
the  model  is  accepted. 


Table  11.  Results  of  24  Tests,  First  Order  Model  of  3  Variables 
Life  predicting  equation: 

y  =  3.56  +  0.0195x,  -  0.0325x„  -  0.0518xo 
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Table  12.  Analysis  of  Variance  of  24  Tests 
First  Order  Model  of  11  Variables 
Life  predicting  equation: 


y  =  3.50  -  0.0105x^  -  0.398x2  -  0.0748x3  +  0.0459x4 

-  O.lllx,  -  0.317x,  +  0.290x,  +  G.C091xo  -  0.0383xo 
jo/  o  9 

+  0.0739x1q  +  0.0231x11 


Source 

Sum  of 
Squares 

Degrees  of 
Freedom 

Mean 

Squares 

F-Ratio 

Due  to  Mean 

0.1223 

1 

0.1223 

Due  to  Variance 

2.3905 

1 

2.3905 

Due  to  Zero  Upcrossings 
Due  to  Level  Upcrossings 

0.0482 

1 

0.0482 

of  Level 

0.0006 

1 

0.0006 

of  e  Level 
y 

0.0838 

1 

0.0838 

Due  to  Duration  of  Excur¬ 
sion  Above  Zero  Level 
Due  to  Duration  of  Excur- 

0.0180 

1 

0.0180 

sion  Above  Level 

0.1010 

1 

0.1010 

Above  e  Level 
y 

0.0006 

1 

0.0006 

Due  to  Band  Width 

Due  to  Average  Amplitude 

0.0238 

1 

0.0238 

Above  ef  Level 

0.0362 

1 

0.0362 

Above  e  Level 

y 

0.0054 

1 

0.0054 

Due  to  Regression 

2.8302 

11 

0.2573 

Residuals 

0.3498 

12 

0.0292 

8.812 

Total 

3.1801 

23 

F-ratio  is  greater  than  the  table  value  of  2.73  with  11  and  12  degrees 
of  freedom  at  95%  significance  level.  So  regression  is  effective,  and 
the  model  is  accepted. 
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Table  14  .  Comparison  of  Percent  Deviations  of  Predicted  Lives 
and  Residual  Sum  of  Squares  for  Five  Models 


Test 

No. 

Actual 
Life  T 

Percent  Deviations  of 

Predicted  Lives 

Eq(l) 

Eq(2) 

Eq(3) 

Eq(4) 

Eq(5) 

25 

49.58 

-3.7 

-3.1 

-10.2 

-2.9 

-9.0 

26 

50.75 

14.9 

11.9 

-11.6 

9.9 

-0.7 

27 

18.41 

-26.4 

-25.7 

-2.0 

-28.9 

-9.6 

28 

52.50 

6.3 

8.1 

7.1 

7.2 

5.5 

29 

22.33 

-15.4 

-17.2 

-15.9 

-18.0 

-10.6 

30 

30.67 

-52.3 

-60.1 

-7.7 

-61.7 

-22.2 

31 

24.00 

2.7 

-5.2 

-5.7 

-6.7 

-5.8 

32 

30.00 

10.4 

8.5 

13.7 

8.4 

6.0 

33 

40.92 

2.1 

19.0 

6.9 

17.4 

3.6 

34 

44.42 

21.5 

18.5 

16.2 

-22.3 

23.2 

35 

38.17 

9.2 

-2.0 

8.5 

5.3 

36 

36.92 

7.3 

2.1 

5.9 

2.4 

37 

33.08 

-2.1 

0.9 

-4.0 

-7.9 

38 

23.33 

13.3 

8.2 

12.5 

4.4 

39 

30.83 

-5.9 

-3.7 

-8.3 

-12.7 

40 

32.17 

-11.5 

-11.4 

-12.1 

-26.0 

41 

40.42 

-1.1 

-7.5 

-0.8 

-17.9 

42 

71.33 

10.5 

14.6 

9.7 

10.2 

43 

19.75 

-27.3 

-2.4 

44 

29.92 

15.4 

19.8 

45 

48.17 

-4.8 

3.1 

46 

44.50 

15.3 

14.2 

47 

21.00 

9.5 

2.2 

48 

64.92 

5.6 

8.8 

Average  Deviations 


Negative  side 
Positive  side 

24.5 

9.7 

14.7 

11.8 

7.8 

8.7 

Residuals 

16.5 

10.4 

13.4 

8.4 

Percent  residual 

sum  of  squares 

of  the  total 

32.6 

24.0 

8.3 

20.5 

11.0 

31 
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Fig.  2  Distribution  of  Residuals  for  the  Best  First  Order  Model,  Eq.  (3) 
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designs  was  conducted  for  24  specimens  for  fatigue  of  materials  in  an  elastic- 
plastic  range  under  random  vibrations.  The  magnitude  levels  of  all  probabi¬ 
listic  parameters  are  completed  from  the  spectral  moments  and  coded.  By 
regressing  the  log  of  fatigue  life  on  the  corresponding  coded  parameters,  the 
first-order  regression  models  are  developed. /m\ 
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